Local response of topological order to an external perturbation 
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We study the behaviour of the Renyi entropies for the toric code subject to a variety of different perturba- 
tions, by means of 2D DMRG and analytical methods. We find that Renyi entropies of different index a display 
derivatives with opposite sign. This splitting happens at a low value of a, as opposed to typical symmetry break- 
ing states, and can be detected on a very small subsystem regardless of the correlation length. This phenomenon 
is due to the presence in the phase of a point with flat entanglement spectrum, zero correlation length, and area 
law for the entanglement entropy. We argue that this kind of splitting is common to all the phases with a certain 
group theoretic structure, including quantum double models, and cluster states. The fact that the size of the 
subsystem does not need to scale with the correlation length, makes it possible for this effect to be accessed 
experimentally. 



The understanding of collective behaviour arising from mi- 
croscopic interactions is one of the central issues in physics. 
Indeed, the most recent research has been indicating that the 
kind of emerging order in extended systems is not so simple 
as was expected. This is the case of quantum phases of mat- 
ter with the so-called topological order. These phases, and 
the phase transitions between them, cannot be characterized 
within the symmetry breaking mechanism, a cornerstone in 
many-body physics. Topological phases are characterized by 
global correlations that are captured by subtle entanglement 
properties. This makes topological order an elusive kind of 
order, and therefore very difficult to detect. Here we show 
that the order in a paradigmatic class of topological states can 
be revealed by measurements on a small portion of the system 
through the analysis of their response to an external perturba- 
tion. To this aim, we employ analytical and numerical meth- 
ods in quantum many-body theory and quantum information. 
Our results can provide a venue to detect topological order 
experimentally. 

The Landau symmetry breaking mechanism classifies dif- 
ferent phases of many-body systems according to the symme- 
try that the system-low-energy-states breaks [1]. This implies 
that some observable exists whose values tell us in what phase 
the system is, the so called order parameter. In this picture, the 
macroscopic order arises from local order. For instance, we 
can distinguish a magnet from a paramagnet by measuring the 
magnetization resulting from the breaking of a local rotation 
symmetry of the atomic spin. 

In the last two decades, a more complex scenario emerged. 



An example of paramount importance in condensed matter 
physics is provided by quantum phase transitions occurring 
in two dimensional electron gas displaying fractional quan- 
tum Hall effect[2] : Different Hall states are characterized by 
different physical properties, despite sharing the same sym- 
metries. Different types of topological orders have been found 
in various physical systems so far, like spin liquids [3], any- 
onic systems [4], topological insulators [5]. Since these states 
do not break any symmetry, they cannot be distinguished 
by measuring a local order parameter. The comprehension 
of topological phases is one of the central issues in mod- 
ern condensed matter[4]. Topological states have properties 
like protected edge modes or degeneracies that are robust un- 
der arbitrary perturbations [4]. They possess anyonic exci- 
tations whose interactions are of topological nature and are 
inherently robust against decoherence. These features have 
made such states very interesting as a platform for quantum 
computers [6]. 

As a result from the interaction between condensed matter 
and quantum information, it has been understood that topo- 
logical order in the ground state of a given system can be 
characterized through entanglement[7-9]. Specifically, a 
long range pattern of entanglement results in the ground state, 
called topological entanglement entropy j[9]. Such a quan- 
tity can be computed through the von Neumann entropy Si 
of the reduced density matrix pa describing the subsystem A 
(being the system A u B in pure state). For a regular region 
A then von Neumann entropy is given by Si = a\dA\ - 7 + g, 
where \dA\ measures the boundary of A. The quantities a, g 
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FIG. 1. Cylinder of infinite length and width L y = 5 used in 2D 
DMRG calculation. (A) Subsystems on which Renyi entropies are 
calculated: A a - one star and A SiP - composition of star and pla- 
quette. Loops l\ and if used to distinguished between topological 
sectors are also depicted. (B) Subsystem Ch i. that contains half of 
the infinite cylinder . 

are non universal and depend on the details of the Hamiltonian 
describing the system, like the correlation length £. The finite 
correction 7, instead, is universal and can serve as partial clas- 
sification of topological phases[10]. We remark that, although 
A is a finite region, Si is not a local order parameter, because 
it is not the expectation value of a local observable. On the 
other hand, the reconstruction of pa by local measurements 
seems to be quite a daunting task because A, though finite, 
scales with £[9]. Larger system sizes can be studied numeri- 
cally by using Renyi entropies through the replica trick[l 1]. 

Much effort has been devoted to analyze the topological or- 
der in the ground state for a fixed value of the control parame- 
ter A. In this paper, we push forward the idea that progress in 
probing topological order can be achieved by analyzing local 
variations of the ground state by external perturbation. 

We study a set of spin- 1/2 localized at the edges of 2D 
square lattice with periodic boundary conditions in presence 
of a perturbation V: 

^ = -En<-EEK + nA) (i) 

s its p iep 

where s and p label the vertices and plaquettes of the lattice 
respectively, while erf , of are Pauli operators of the spin liv- 
ing at the edge i. For V^(A) = the Hamiltonian above is the 
celebrated toric code, a paradigmatic model for topological 
order[12]. For the analysis below, we remark that in this case, 
the ground state of this model features £ = 0. We consider 
different V(X) (see Table I) where A stands for {Ai, . . . , A^} 
that are the parameters controlling the perturbation. The per- 
turbation in (1) such that the correlation length is increasing 
with Ai until divergence at the critical point A c . 

For A < A c these systems are topologically ordered, while 
for A > A c they are trivial paramagnets. In both phases there 
is no local order parameter. This model belongs to a class of 
so-called quantum double models which correspond to those 
phases whose low energy theory is a lattice gauge theory [12]. 




FIG. 2. The splitting phenomenon. The figure displays the splitting 
with opposite slopes between the small and large a Renyi entropies. 
We see the splitting occurring around a ^ 0.6. The Renyi entropies 
are calculated for the partition A StP of Fig. 1 A for the ground state of 
H Htc ■ \\ . 

In this article, we demonstrate how the topological phase can 
be distinguished from the paramagnet of (1) through the in- 
spection of d\S a , where S a = (1 - a) -1 logTrp^ are the 
Renyi entropies. We shall see that d\S a present a peculiar 
behavior even for a small subsystem A as £ varies as an effect 
of a perturbation. We remind that So = log R, where R is the 
Schmidt rank of the state, i.e. the number of nonzero eigenval- 
ues of pa, and that Si is the von Neumann entropy measuring 
the entanglement entropy for the subsystem A. We observe 
that, generically, R and thus So increase with £ because more 
degrees of freedom get entangled by increasing the room for 
correlations. On the other hand, we find that the behaviour of 
Si is drastically different between topological phases and the 
topologically trivial phases presented here. 

Our approach elaborates on the quantum information no- 
tion of differential local convertibility (DLC), which states 
that two proximal bipartite pure states are locally convertible 
if it is possible to transform one into the other by resorting 
only on local quantum operations on A and B separately, plus 
classical communication and the aid of an ancillary entangled 
system[13]. The necessary and sufficient condition for that 
is d\S a < for all a > [14]. In Ref.[15] such notion has 
been applied to the order-disorder quantum phase transition 
occurred in a one dimensional spin system. 

For each V(X) we compute the ground state wavefunction 
\ip(X)) and its reduced density matrix pa(X). For some V^(A) 
we can apply exact analytical approach; for the generic per- 
turbation V XZ (X) = T,i(X z o~i + X x af) we resort to numerical 
analysis. 

The numerical method employed here is an infinite DMRG 
algorithm [29] in two dimensions, detailed in [26]. The 
method provides Matrix Product State (MPS) representation 
of a complete set of ground states on a cylinder of infinite 
length and finite width L y (Fig.l) for a given Hamiltonian 
that realizes topological order. As argued in [28], each ground 
state has a well-defined flux threading through the cylinder. 
The flux is measured by (in general) dressed Wilson loop op- 
erators that enclose the cylinder in the vertical direction. 

In the case of fixed point toric code (Eq.l with V = 0), 
these loops are given by if and if (Fig.l A). Four topological 
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TABLE I. Local Convertibility in Topological Phase. DLC, i.e no 
splitting of the Renyi entropies, only occurs if the perturbation is fine 
tuned in order to keep the system with £ = 0. The left column shows 
the type of perturbation studied. The first column details whether the 
considered model is Gauge Invariant. The second column indicates 
whether DLC occurs. For certain perturbations the ground state of 
the system is accessible exactly (third column). The last column pro- 
vides the information on f . The subsystems we refer here have a non 
trivial bulk[21] . 



sectors are then distinguished by (if), (if) = ±1- Once the 
perturbation is present, Wilson loops may change, but as long 
as the perturbation is small, (Zf) and (if) can still be used to 
identify topological sectors because (Zf ), (if) - ±1. 

Simulations are carried out with cylinders of width up to 
L y = 5 for sj\ 2 x + X 2 Z < 0.05 and < A < 0.7 as shown in 
Fig. 2 and Fig. 3 respectively. In the topological phase, the 
outcome of each simulation is four quasi-degenerate ground 
states, from which the one with (Zf),(/f) - + 1 i s chosen 
for further investigation. This is done to ensure that finite 
size effects have the least possible impact on results. In the 
limit L y -* oo all four ground states become locally indistin- 
guishable. The results are converged in bond dimension of 
MPS which acts as a refinement parameter. A reduced density 
matrix of a half-infinite cylinder Ch.i (Fig. IB) is calculated 
throughout the simulation. The bond dimension is increased 
until convergence of its spectrum is reached. 

In Fig. 2, we can see the behaviour of the S a Renyi entropies 
as we span the parameter space X X ,X Z for the perturbation 
V xz . We see clearly that in the topologically ordered phase a 
splitting of S a 's occur: d\S a £ at a given value of a = a$; 
we found o.q - 0.6 (see caption of 2). We call this phe- 
nomenon low-a splitting in the rest of this article. In the para- 
magnetic phase all the Renyi entropies are monotones with 
A. This behaviour is generically independent of the size and 
shape of the subsystem A, as long as A contains some bulk 
[21]. Below we provide an explanation of the phenomenon. 
The topologically ordered phase we consider is characterized 
by the presence of a state (at A = 0) with £ = and a flat 
entanglement spectrum (and an area law) [10]. The flat entan- 
glement spectrum implies that small perturbations results in 
decreasing S a for a > a being a < 1, because the distribu- 
tion becomes less flat in the most represented eigenvalues in 
the entanglement spectrum. In contrast, So must increase with 
£ as an effect of the perturbation (new degrees of freedom are 
involved in the entanglement spectrum). So the low-a split- 
ting results from the insertion of a finite £ in the state evolving 
from a state with flat spectrum and zero £. We also observe 
that such property is shared with the so called G-states which 




0.1 0.2 0.3 \ 0.4 0.5 0.6 0.7 



FIG. 3. Renyi entropies as a function of A for the ground state of 
H = Htc + Vxz with X x = X and X z = A/2. Here, L y = 5. The re- 
duced system A consists of A s and Ch.i. in panels A and B, respec- 
tively. As A increases, the correlation length increases. The Schmidt 
rank R and the low a < ao-Renyi entropies increase as well. The 
value of ao are 0.4, 0.6 in panel A and B, respectively. Nevertheless, 
the total entanglement Si and all the higher Renyi entropies are de- 
creasing with £. Notice the spike in panel B marking the quantum 
phase transition to the paramagnetic phase at A c ~ 0.35. 

include all the topologically ordered quantum double models, 
and states like the cluster states [27], and therefore our find- 
ings apply to this class of models as well[10] (see [31] for a 
discussion of the cluster phase diagram). For the correspond- 
ing ground states the splitting is thus also expected to occur 
at ao < 1 no matter how we perturb in the plane X X ,X Z . Note 
again that in the paramagnetic phase all the d\S a have the 
same sign and no splitting ever occurs, which is easily un- 
derstood from the presence (at very large A) of a completely 
factorized state, see Fig. 3. 

Now we argue how the specific values at which the split- 
ting phenomenon occurs effectively distinguish the topolog- 
ically ordered state from a topologically- trivial ordered state 
(like a ferromagnet). Indeed, as discussed above, the latter 
states have typically Si increasing with £ (S2 is also increas- 
ing as well as long as £ is not too large). Accordingly, the 
splitting phenomenon is shown to occur at higher values of a 
for phases with broken symmetries [15]. We also observe that 
it may or may not occur for these phases, depending on local 
details of the system (the anisotropy in the case of the study 
in Ref.[15]). 

Summarizing: We can distinguish between the topological 
phase and the paramagnet of (1). Furthermore, we can distin- 
guish between the topological phase and a symmetry break- 
ing phase because the splitting phenomenologically occurs at 
a > 1 in the latter. 

To corroborate our findings, we resort to exact analysis 
for suitable perturbations V(X)'s. We consider two cases 
i) Vh = X z Y,ieh a i> corresponding to placing the external 
field oc cr z only along the horizontal links of the lattice; ii) 
V(X) = 2 s e A ^" s<Ti leading to the Castelnovo-Chamon 
model[20]. Since these perturbations commute with the pla- 
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quette operators of Eq. 1 , the ground state of these models can 
be written as the superposition of loop states \g) with ampli- 
tudes a(g). A loop state \g) is obtained from the completely 
polarized state in the z direction, by flipping down all the spins 
intersected by a loop drawn on the lattice. The corresponding 
loop operators g form a group G called the gauge group of 
these theories. 

In the case i) the star operators Ylies a i interact only along 
the rows of the lattice. The model maps onto the product of 
arrays of Ising chains by the duality A s -* rf, erf -* t^t£ +1 : 
Htc + V h m. H f/ = efciC-A^ t*t* +1 - Z^rp [17, 18]. 
The relevant correlators in the variables a can be obtained 
through the correlators in the dual variables r that can be ac- 
cessed exactly[19]. In the following, we sketch a proof that 
the splitting phenomenon does occur in this model (see [2 1 ] 
for details). We consider the star A s = {i%,i2,isi h} as sub- 
system A (see Fig. 1); pa, is block diagonal with 2x2 blocks 
labeled by \i1i2i3i4) and A s ]ii«2*3*4)- It results that pa has 
maximum rank unless a(g) = a(gi)a(g2), implying there is 
a zero eigenvalue in each block. In the dual picture this is 
equivalent to require (riTj) = (t,;)(t.,). Such condition holds 
at A = only, and therefore R increases at A * 0. The factor- 
ization of the amplitudes also proves that both the a = 1,2- 
Renyi entropies decrease at small A[21]. 

The case ii) is important to test the argument of the interplay 
between splitting and increasing of correlation length. This 
argument implies that a perturbation for which £(A) = const 
does not lead to a splitting in the Renyi entropies. The model 
of Castelnovo-Chamon features exactly this, since spin-spin 
correlation functions (afcr*) are vanishing for every value of 
A. The exact ground state is made of loops with amplitudes 
a(g) = e -V2E*«<(9) i where a 7 ({g) = {g\a z t \g). The topo- 
logical phase is(A < A c « 0.44). A lengthy calculation leads 
toS a (p A ) = (l-a)- 1 logZ- Q (A)Z geG e- Ai ^«- 1 (A, 3 ), 
where Z = Y g e~ XL * and w(\,g) := Z h , GAMGB e~ XL ^, 
L m is the length of the loop m of the gauge group G; here 
Ga and Gb are the gauge groups of the subsystems A and 
B respectively. The analysis of small and large A expansions 
reveals that d\S a (X) < 0[21]. As a particular case, Sq is con- 
stant for every value of A. Accordingly, for this fine tuned per- 
turbation all Renyi entropies decrease and therefore no split- 
ting is observed. This is consistent with the fact that also in 
this model the amplitudes a(g) factorize as discussed in i) . 

In this paper, we have shown that topologically ordered 
states in the paradigmatic class of quantum double models 
have a peculiar way of responding to perturbations: The 
Renyi entropies of a subsystem display opposite slopes around 
a ^ 0.6. The subsystem size can be small independently of 
£. In the paramagnetic phases, in contrast, no such a phe- 
nomenon is observed. The low-a splitting is shown to be a 
robust characteristic of the topological phase, see Fig. 2. In a 
quantum information context, we can say that, while param- 
agnets are always differential locally convertible, the class of 
topologically ordered ground states we studied here are not. 

We close this paper with some speculations and future di- 



rections. Comparing our findings with the study in [15], we 
see that the kind of Renyi entropies splitting that may oc- 
cur in symmetry-breaking-phases and the one characterizing 
topological states are phenomenologically very different. We 
conjecture that in topologically-trivial ordered phases the low- 
a-splitting can only be accidental, meaning that every sym- 
metry breaking state can be continuously "deformed" with- 
out observing it. Therefore, by measuring £ together with 
the von-Neumann entropy Si for a small subsystem (it can 
be shown this is sufficient to show whether low-a splitting is 
present), one could tell apart topological states of the class 
studied here from symmetry breaking states. This should be 
contrasted with the standard approach based on topological 
entropy where the size of the subsystem must be much larger 
than £. 

Our results may provide the basis for the experimental de- 
tection of exotic phases like topologically ordered or cluster 
states. In this endeavour, it would be desirable to investigate 
these findings in presence of disorder. 

It would be also important to investigate the behaviour of 
Renyi entropies under perturbations for more general topolog- 
ical states without a flat entanglement spectrum like fractional 
quantum Hall liquids and chiral spin liquids[ll, 24], and for 
symmetry protected topological states like one dimensional 
systems or topological insulators[25]. A relevant question is 
whether the splitting behaviour can shed light on the finite 
temperature resilience of topological order, and on the very 
existence of anyonic excitations. 

Acknowledgments. — We thank J.I. Cirac, R. Fazio, G. Vi- 
dal, X.G. Wen for discussions. This work was supported 
in part by the National Basic Research Program of China 
Grant 201 1CBA00300, 2011CBA00301 the National Natu- 
ral Science Foundation of China Grant 61073174, 61033001, 
61061130540. Research at Perimeter Institute for Theoreti- 
cal Physics is supported in part by the Government of Canada 
through NSERC and by the Province of Ontario through MRI. 



[1] N. Goldenfeld, Lectures on phase transitions and the renormal- 
ization group, (Addison Wesley, New York, 1992). 

[2] H. L. Stormer, Rev. Mod. Phys. 71, S298S305 (1999). 

[3] S. Yan, D. A. Huse, and S. R. White, Science, 332, 1 173 (201 1). 

[4] X.-G. Wen, Quantum Field Theory of Many-body Systems (Ox- 
ford University Press, 2004). 

[5] M. Z. Hasan and C. L. Kane, Rev. Mod. Phys. 82, 3045 (2010). 

[6] M. H. Freedman, A. Kitaev, and Z. Wang, Commun. Math. 
Phys. 227, 587 (2002); C. Nayak, S. H. Simon, A. Stern, M. 
Freedman, and S. Das Sarma, Rev. Mod. Phys. 80, 1083 (2008). 

[7] L. Amico et al., Rev. Mod. Phys. 80, 517 (2008). 

[8] J. Eisert, M. Cramer, and M. B. Plenio, Rev. Mod. Phys. 82, 
277 (2010). 

[9] A. Hamma, R. Ionicioiu, and P. Zanardi, Phys. Lett. A 337, 22 
(2005); A. Hamma, R. Ionicioiu, and P. Zanardi, Phys. Rev. A 
71, 022315 (2005); A. Kitaev and J. PreskiU, Phys. Rev. Lett. 
96, 110404 (2006); M. Levin and X.-G. Wen, Phys. Rev. Lett. 
96, 110405 (2006). 



5 



[10] S. T. Flammia, A. Hamma, T. L. Hughes, and X.-G. Wen, Phys. 

Rev. Lett. 103, 261601 (2009). 
[11] S. V. Isakov, M. B. Hastings, and R. G. Melko, Nature Phys. 7, 

772(2011). 

[12] A. Y. Kitaev, Ann. Phys. (N. Y.) 303, 2 (2003). 

[13] D. Jonathan and M. Plenio, Phys. Rev. Lett. 83, 3566 (1999). 

[14] S. Turgut, J. Phys. A: Math. Theor. 40, 12185 (2007); M. 

Klimesh, arXiv:0709.3680 (2007). 
[15] J Cui, M Gu, LC Kwek, MF Santos, H Fan, V Vedral, Nat 

Commun 3, 812 (2012). 
[16] G. Vidal, J. I. Latorre, E. Rico, A. Kitaev, Phys. Rev. Lett.90 

227902 (2003). 

[17] J. Yu, S.-P. Kou, and X.-G. Wen, Europhys. Lett. 84, 17004 
(2008). 

[18] G. B. Halasz and A. Hamma, arXiv: 1206.4223. 
[19] E. Barouch and B. M. McCoy, Phys. Rev. A 3, 786 (1971). 
[20] C. Castelnovo and C. Chamon, Phys. Rev. B 77, 054433 (2008). 
[21] Some additional care is to be used when considering a subsys- 
tem A which is all boundary Ga = {Ia}. A. Hamma, et al. in 



preparation. 

[22] D. Abanin, and E. Dernier, Phys. Rev. Lett. 109, 020504 (2012). 
[23] X.-G. Wen, Adv. Phys. 44, 405 (1995). 

[24] S. Papanikolaou, K. S. Raman, and E. Fradkin, Phys. Rev. B 76, 

224421 (2007). 
[25] J. Cui, L. Amico, M. Gu, and V. Vedral, in preparation. 
[26] L. Cincio, G. Vidal, in preparation. 

[27] A. Hamma, R. Ionicioiu, P. Zanardi, Phys. Rev. A 72, 012324 
(2005) 

[28] L. Cincio, G. Vidal, arXiv: 1208.2623. 

[29] S. R. White, Phys. Rev. Lett. 69, 2863 (1992); I. P. McCulloch, 
arXiv:0804.2509; G. M. Crosswhite, Phys. Rev. B 78, 035116 
(2008). 

[30] A. Hamma, W. Zhang, S. Haas, and D. A. Lidar, Phys. Rev. B 

77,155111 (2008). 
[31] H. Kalis, D. Klagges, R. Orus, K.P. Schmidt, Phys. Rev. A 86, 

022317 (2012) 

[32] P. Horodecki and A. Ekert, Phys. Rev. Lett. 89, 127902 (2002); 

Phys. Rev. Lett. 108, 110503 (2012) 
[33] We expect the S2 test to be reliable only for moderate £. 



